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Abstract — Monotonicity criteria are established for the gener- 
alized Marcum g-function, Q M (a, (3), the standard Nuttall g- 
function, Q M N (a, (3), and the normalized Nuttall g-function, 
Qm,n{ch, (3), with respect to their real order indices M,N. 
Besides, closed-form expressions are derived for the computation 
of the standard and normalized Nuttall g-functions for the case 
when M, N are odd multiples of 0.5 and M > N. By exploiting 
these results, novel upper and lower bounds for Q M N (a, f3) 
and Qm,n{ol, f3) are proposed. Furthermore, specific tight 
upper and lower bounds for Q M (a, (3), previously reported in 
the literature, are extended for real values of M . The offered 
theoretical results can be efficiently applied in the study of 
digital communications over fading channels, in the information- 
theoretic analysis of multiple-input multiple-output systems and 
in the description of stochastic processes in probability theory, 
among others. 

Index Terms — Closed-form expressions, generalized Marcum 
g-function, lower and upper bounds, monotonicity, normalized 
Nuttall g-function, standard Nuttall g-function. 



I. Introduction 
A. The Nuttall Q-Functions 

N extended version of the (standard) Marcum Q- 



. function, Q(a,f3) = J 



2 I Q (ax)dx, where 
a, f3 > 0, originally appeared in (U Appendix, eq. (16)], 
defines the standarqj Nuttall (2-function |2, eq. (86)], given 
by the integral representation 



Qm,n{(x,P) 



x M e 



' lN{ax)dx 



(1) 







where the order indices are generally reals with values M > 
and N > — 1, 7jy is the iVth order modified Bessel function 
of the first kind |[3] eq. (9.6.3)] and a, (3 are real parameters 
with a > 0, (3 > 0. It is worth mentioning here, that the 
negative values of N, defined above, have not been of interest 
in any practical applications so far. However, the extension 
of the Nuttall g-f unc tion to negative values of N has been 
introduced here in order to facilitate more effectively the 
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relation of this function to the more common generalized 
Marcum g-function, as will be shown in the sequel. An 
alternative version of Q M N (a, (3) is the normalized Nuttall 
Q-function, Qm,n(&, 0), which constitutes a normalization of 
the former with respect to the parameter a, defined simply by 
the relation 



IM.N 



(a, 13) 



(2) 



Typical applications involving the standard and normalized 
Nuttall Q-functions include: (a) the error probability per- 
formance of noncoherent digital communication over Nak- 
agami fading channels with interference ID, (b) the outage 
probability of wireless communication systems where the 
Nakagami/Rician faded desired signals are subject to indepen- 
dent and identically distributed (i.i.d.) Rician/Nakagami faded 
interferers, respectively, under the assumptions of minimum 
interference and signal power constraints [4|-[7|, (c) the per- 
formance analysis and capacity statistics of uncoded multiple- 
input multiple-output (MIMO) systems operating over Rician 
fading channels [8 ]-[ 10], and (d) the extraction of the required 
log-likelihood ratio for the decoding of differential phase-shift 
keying (DPSK) signals employing turbo or low-density parity- 
check (LDPC) codes ifTTl . 

Since both types of the Nuttall Q-fxmc\ion are not consid- 
ered to be tabulated functions, their computation involved in 
the aforementioned applications was handled considering the 
two distinct cases of M + N being either odd or even, in 
order to express them in terms of more common functions. The 
possibility of doing such when M + N is odd was suggested 
in 0, requiring particular combination of the two recursive 
relations |j2] eqs. (87), (88)]. However, the explicit solution 
was derived only in fl4] eq. (13)] entirely in terms of the 
Marcum Q-function and a finite weighted sum of modified 
Bessel functions of the first kind. Having all the above in 
mind, along with the fact that the calculation of Q(a,f3) 
itself requires numerical integration, the issue of the efficient 
computation of ([TJ and (O still remains open. 



B. The Generalized Marcum Q-Function 

The generalized Marcum g-function [12| of positive real 
order M, is defined by the integral lfl3l eq. (1)] 



Gm(",/3) = 



vAT-l 



M - g +° 

x e 2 



" lM-i{oix)dx (3) 
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where a, (3 are non-negative real parameters^. For M = 1, 
it reduces to the popular standard (or first-order) Marcum Q- 
function, Q 1 (a, /3) (or Q(a, f3)), while for general M it is 
related to the normalized Nuttall Q-function according to [14 
eq. (4.105)] 

2jf(a^) = SM,M-i(a,/3) ) a>{). (4) 

An identical function to the generalized Marcum Q is the 
probability of detection^] [1 1 eq. (49)], which has a long history 
in radar communications and particularly in the study of target 
detection by pulsed radar with single or multiple observations 
ID, ASHES. Additionally, Q M (a, f3) is strongly associated 
with: (a) the error probability performance of noncoherent and 
differentially coherent modulations over generalized fading 
channels fl4l . Ifl9l - ll23l . (b) the signal energy detection of 
a primary user over a multipath channel ll24l . 11251 . and finally 
(c) the information-theoretic study of MLMO systems |26|. 
Aside from these applications, the generalized Marcum Q- 
function presents a variety of interesting probabilistic interpre- 
tations. Most indicatively, for integer M, it is the complemen- 
tary cumulative distribution function (CCDF) of a noncentral 
chi-square (\ 2 ) random variable with 2M degrees of freedom 
(DOF) 11271 eq. (2.45)]. This relationship was extended in 
1 28 1 to work for the case of odd DOF as well, through a 
generalization of the noncentral \ 2 CCDF. Similar relations 
can be found in the literature involving the generalized Rician 
||29l (2.1-145)], the generalized Rayleigh |f30l pp. 1] (for 
a = 0) and the bivariate Rayleigh [31, Appendix A], 11321 
(for M = 1) CCDF's. Finally, in a recent work IT3311 . a new 
association has been derived between the generalized Marcum 
2-function and a probabilistic comparison of two independent 
Poisson random variables. 

More than thirty algorithms have been proposed in the 
literature for the numerical computation of the standard and 
generalized Marcum (^-functions, among them power series 
expansions [34|-[36|, approximations and asymptotic expres- 
sions lT37ll - ll40ll . and Neumann series expansions BH -I43I. 
However, the above representations may not always provide 
sufficient information about the relative position of the ap- 
proximated value with respect to the exact one, which in some 
applications is highly desired. In 11441 . the generalized Marcum 
2-function of integer order M has been expressed as a single 
integral with finite limits, which is computationally more 
desirable relatively to other methods suggested previously. 
Nevertheless, the integral cannot be computed analytically 
and appropriate numerical integration techniques have to be 
applied, thereby introducing an approximation error in its 
computation. In 11451 . an exact representation for Q M (a,/3), 
when M is an odd multiple of 0.5, was given as a finite sum 
of tabulated functions, assuming that f3 2 > a 2 + 2M. This 
result was recently enhanced in [46 1 to a single expression 
that remains accurate over all ranges of the parameters a, j3, 
while in l47l the same expression was bounded by particular 

2 For a = the right hand side of {5) can be easily shown to satisfy the 
limiting value of 1141 eq. (4.71)], reproduced in i30\ . 

3 For M incoherently integrated signals, the two functions are simply related 
by Q u {a,P) = PmOm' If)' as induced ^ El eq. (7)]. 



utilization of novel Gaussian (^-function inequalities. Finally, 
in (48l , an equivalent expression to ll46l eq. (11)] was derived, 
adopting a completely different (analytical) approach from the 
latter. 

Close inspection of the issues mentioned above, render the 
existence of upper and lower bounds a matter of essential 
importance in the computation of (0. Several types of bounds 
for the standard Il49l-ll5ll and generalized EJ, (46), (48l , 
[52 1 — [ 54- 1 Marcum 2-functions nave been suggested so far. 
However, all the aforementioned works consider just integer 
values of M, which is generally true when this parameter 
represents the number of independent samples of a square- 
law detector output. Nevertheless, in many applications, this 
requirement does not hold. According to [14, Sec. 4.4.2], it 
would be desirable to obtain alternative representations for 
Q M (a,/3) regardless of whether M is integer or not. For 
instance, in l20ll . the fading parameter of the Nakagami- 
m distribution is restricted to integer values in the lack of 
a closed-form expression for the generalized Marcum Q- 
function of real order. Additionally, in (25), (55)-|l57), the 
order M of the generalized Marcum (^-function, involved in 
the energy detection in various radiometer and cognitive radio 
applications, is expressed as the product of the integration 
time and the receiver bandwidth, thus implying that in general 
M is a non-integer quantity. Furthermore, a probabilistic 
interpretation of Q M _ /1 (a, /3), where M € and fj, = 0.5, 
is given in (58), where it is related to main probabilistic 
characteristics of 2(M — /i) random variables, while in [59 1- 
(6T1 . noncentral \ 2 random variables with fractional DOF are 
studied. Finally, in [62 1, the integrand of (O has been proved 
to be a probability density function (PDF) for a > and 
M > 0, a result that also has been utilized in (48). 

C. Contribution 

As described in Subsection II-AI a closed-form expression 
for the computation of the standard and normalized Nuttall 
(^-functions is available in the literature only for the case of 
odd M + N, with the additional restriction of integers M, N. 
In Subsection III-AI we derive a novel closed-form expression 
for the computation of Qm Ni a i P) an d Qm,n(&, (3) when 
M, N are odd multiples of 0.5 and M > N, being valid for 
all ranges of the parameters a, (3. 

Besides, in Subsection IH-BI we proceed with the estab- 
lishment of appropriate monotonicity criteria, revealing the 
behavior of both functions with the sum M + N. Specifically, 
we demonstrate that the standard Nuttall Q-fimctiori is strictly 
increasing with respect to M+N when M > 7Y+1, under the 
constraints of a > 1 and j3 > 0. For the normalized Nuttall Q- 
function, a similar monotonicity statement is proved without 
the necessity of reducing the range of a. 

An alternative approach, sufficient enough to facilitate the 
problem of evaluating the Nuttall Q-functions, is the derivation 
of tight bounds. Nevertheless, to the best of the authors' 

4 Throughout the manuscript, we adopt IN and Wo notations for the 
representation of the positive and the non-negative integer set, respectively. 
In the same way, R+ includes the positive and Rq*" the non-negative reals. 
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knowledge, such bounds have not been reported in the litera- 
ture so far. Subsection III-B I is completed with the exploitation 
of the previous results in order to derive novel upper and lower 
bounds for Q M N (a, j3) and Qm,n(&, 0) when M > N + 1 
and f3 > 0, with the extra requirement of a > 1 for the former. 

Additionally, in Subsection II-BI the need for computing the 
generalized Marcum g-function, Q M (a,/3), of real order M 
was highlighted, since it is a case of frequent occurrence in 
various applications. However, a thorough literature search 
for studies concerning arbitrary values of M, revealed only 
[46| for the closed-form computation of Q M (a,f3) of half- 
odd integer order, and the accepted paper [63], where bounds 
for Q M (a,f3) were introduced for the case when M is not 
necessarily an integer. These considerations motivated us to 
generalize the scope of M in ll46l eq. (16)] 

Qm-o.&{<*,0) < Qm{*,P) < Gjif+ojs(a.j8). M e IN (5) 

as described in Section [TTTJ by providing a monotonicity 
formalization for the generalized Marcum g-function, namely 
that Q M (a,/3) is strictly increasing with respect to its order 
M > for a > and /? > 0. This interesting statement was 
also recently presented in (62), using a different approach. As 
a consequence, novel upper and lower bounds for Q M (a,(3) 
of positive real order are derived. We finalize the paper with 
some concluding remarks, given in Section [TV] 



II. Monotonicity of the Nuttall ^-Functions 

A. Novel Closed-Form Representations 

So far, closed-form expression for either type of the Nuttall 
g-function is not available in the literature. In this section, 
we derive such a representation for the case when M, N are 
odd multiples of 0.5 and M > N, through the theorem and 
corollary established below. Before proceeding further with the 
corresponding proofs, some definitions of essential functions 
and notations used, would be very convenient. 

Hereafter, T, 7 and r(-,) will denote the Euler gamma 
[3] eq. (6.1.1)], the lower incomplete gamma [3, eq. (6.5.2)] 
and the upper incomplete gamma J3] eq. (6.5.3)] functions, 
respectively, defined by the integrals 

p OO f-X 

T(z) = t^e^dt, j(z,x) = f-^dt 
Jo Jo 
r(z,x)=T(z)-j{z,x), zeR + ,xeR. 

Notations n\, (to)„ and ('") imply the factorial |3] eq. (6.1.6)], 
the rising factorial (Pochhammer's symbol) 13] eq. (6.1.22)] 
and the binomial coefficient JJJ eq. (24.1.1 C)], respectively, 
defined by n! = rifc=i ^ f° r n, G ff; = 1 for n — 0, (m) n = 
( 7 + "~V ! for to e IN, n e Mo and ( m ) = ,. m! ., for 

(m— 1)! ' u \nJ n\(rn — n)\ 

m,n £ f^o, to > n. Finally, sgn(z) = z/\z\ for z ^ 0; = 
for z = 0, stands for the signum function. 

Theorem 1 (Closed-form for the standard Nuttall Q): 
The standard Nuttall 2-fu nc u° n > GmaK^j/O' when 
m = M + 0.5 G M, n = N + 0.5 £ IN and M > N, can 
be evaluated for a > 0, (3 > by the following closed-form 



expression: 



QM,N(a,/3) 



(-l)™(2a)- ,l +2 



(n-k) n -t{2a) k rrk 

X > ^mA^P) 



E 

k=0 



where the term 1^ n (a, (3) is given by 

m—n+k 



1=0 



in - 11 +k\ 2 t^ am _ ll + , , 



(sgn(/3 - a)) ^ 7 1 



2 ' 



l + l 



(6) 



Proof: Given that n = ^ + 0.5 £ IN, the modified Bessel 
function of the first kind, 1^, can be expressed by the finite 
sum [64 eq. (8.467)], which after some manipulations can be 
written as 



In{z) 



(-l)"(2z)-«+3 ^ i (n-fc)„_ 1 (2z) 



E 

fe=0 



k\ 



(l-{-l) k e 2z ), n = N+- Gl.zeR. 



(7) 

Therefore, using ([1} and (0, the standard Nuttall Q-function 
satisfies 

(-l)"(2a)-"+^ ^ (n - fc)„-i(2a) fc 



k=0 

^m-n+feg — dx 



k\ 



m — n+k 



2 dx 



(8) 



-(-ir / » 

h 

The calculation of the integral difference in (HJ can be effec 
tively facilitated by the following definition 



(9) 

where L = m — n + k. Since we examine the case when 
M > iV or equivalently m > n, it follows that in the above 
expression the exponent L is a non-negative integer. Therefore, 
using [65, eq. (1.3.3.18)], © obtains the form 

L 



T k L {a,P)=Y, 



1=0 



,L-l 



{-l) L - 1 / x l e~^dx 



- * 2 dx 



(10) 



-(-1)* / art 

The two integrals involved in ( fTOb can be considered as special 
cases of the more general one 



x e 2 dx, del 
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which for the case of non-negative values of b can be calcu- 
lated from 1641 eq. (3.381.3)] as 

'Z + l 

2~ 

eqs. (3.381.1), (3.381.4)] 



x b 



i— i 

2— r 



6> 



(11) 



while for negative values of b, 
can be combined to yield 



6 < 0. 
(12) 

Therefore, a single expression for the integral X\ for any real 
value of b can be derived, by merging ( fTTT i and ( fT2l with the 
help of [64 eq. (8.356.3)], in order to satisfy 



U =2— 



r 



z + i 



- [sgn(6) 



7 



Z + l 6 



2 \ 1 



Thus, (O is equivalent to 



Z*(a,/?)= £ 



„i-/ 



[(-l) L -^ +Q -(-l) fe I^_ Q ] 



= (-i) fc+1 E(t) 2 ^ a 



. z 

\L-i-fc-l 



Z + l 



-[sgn(/3-a)] i+1 7 



Z + l (^ + a) ; 



2 ' 2 
Z + l (/3-a) 2 



2 ' 2 

which, after the substitution L = m — n + k, yields ©, thus 
completing the proof. ■ 

Corollary 1 (Closed-form for the normalized Nuttall Q): 
The normalized Nuttall (^-function, Qm,n(&, P), when 
to = M + 0.5 G IN, n = N + 0.5 e IN and M > AT, can 
be evaluated for a > 0, (3 > by the following closed-form 
expression: 

(-l) n 2~ n+ i 



IM,N 



(a,/3) 



\fnot 

n-l 

£ 

fc=0 



2n-l 

(n-fc)„_!(2a) fc 



kl 



where the term (a, /3) is given by ©. 

Proof: The proof follows immediately from (|2]i and 
Theorem Q] ■ 

B. Lower one/ Upper Bounds 

In this section, novel lower and upper bounds for the 
normalized and standard Nuttall Q-fimctions are proposed. 
Lemma 1: The function Q s (r,x), defined by 



Gs(r,x) 



a_ T(r + s,x) 

r(r) : 



r, x € H 4 



(13) 



is strictly increasing with respect to r for all s E RJ. 

Proof: By multiplying both the numerator and denomina- 
tor of ( fT3l > by the upper incomplete gamma function, T(r, x), 
we obtain 

JVy -|- g 

Gs(r,x) = — J go(r,z) 



where from ( [TBI one can observe that the term C*o (r, x) is the 
complement of the regularized lower incomplete gamma func- 
tion P(r, x) with respect to unity, defined in [3 , eq. (6.5.1)] by 
P(r, x) = ^ or a ll r > and x e R. Fortunately, P(r, x) 

for r, x > is equal to the cumulative distribution function 
(CDF) of the standard gamma distribution Gamma(r, 1), which 
is strictly decreasing with respect to the shape parameter r. 
Additionally, this important result has also been proved ana- 
lytically in [66, eq. (59)], thus implying that Go{r, x) is strictly 
increasing with respect to r > for all x > 0. Furthermore, 
in ll67l. it has been demonstrated that the function 



K(p,q,x) 



T(p, x) 



p > q > 0,x > (14) 



T(q,x) 

is increasing with respect to q. By substituting p = r + s and 
q = r into ( TT4"1 >. we realize that the ratio T(r + s,x)/T(r,x) 
is increasing with respect to r for s > 0, while it remains 
constant for the trivial case of s = 0. Therefore, it increases 
with r > for all x > 0, s > 0, and the proof is complete. ■ 

The outcome of Lemma Q] will be utilized for the estab- 
lishment of the next theorem, concerning the monotonicity 
property of the normalized Nuttall (^-function. 

Theorem 2 (Monotonicity of the normalized Nuttall Q): 
The normalized Nuttall 2-f unc ti° n ' Qm,n{oi, 0), where 
M > 0, N > — 1 and a,/3 > 0, is strictly increasing with 
respect to the sum M + N, under the requirement of constant 
difference M - N > 1. 

Proof: Combining (Q~|), (0 and using the series represen- 
tation of the modified Bessel function of the first kind in terms 
of the gamma function ll64l eq. (8.445)], we obtain 

Qm,n(<x>P) = e" 



^ k\T(k + N + l)2 2fe " 

oo 



■N 



X 



2k+M+N, 



' dx 



(15) 



where we have interchanged the order of integration and sum- 
mation, since all integrand quantities of the normalized Nuttall 
2-function are Riemann integrable on [j3, oo). Additionally, 
the integral in ( fl3T > is the case of ( fTTT i. thus yielding 



x 2k+M+N e - 



dx 



_ o^ 1 



x r [k + 



M + N + l 



2 



Therefore, ( fT~5b reads 



Qm,n{o-,P) = e c 2 



■2k 



r k 



M+N + l P J 



k=0 



2k+ i^M±l kl T(k + N + l) ■ 

(16) 

Introducing the variables v = M + N and c = M — N and 
taking the partial derivative of both sides of ( [TBI with respect 
to v, we can easily obtain 



d_ 

dv 



-y— 

fe=0 z 

d 



2 k 



-k\ 



du(v) 



G c-i u{v) 



(17) 
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where the function u(v) = k + 1 + 2^ has been employed 
for notational convenience. We note here that, applying the 
Weierstrass M-test [68], the series in ( TTTb can be proved to 
converge uniformly, thus enabling one to interchange the order 
of differentiation and summation. Hence, recalling Lemma Q] 
and the requirement of a > 0, that follows from the definition 
of the normalized Nuttall (^-function, the proof is complete. 



In Figs. U |ta)| and Li lEb)] the normalized Nuttall Q-function has 
been plotted versus the sum M + N for several values of 
a, f3, considering M — N = 1 and M — N = 2, respectively. 
However, we note here that Theorem [2] implies non-integer 
differences M — N as well. 

For the interpretation of the next proposition we define the 
pair of half-integer rounding operators L^Jo.5 an d f^lo.s that 
map a real x to its nearest left and right half-odd integer, 
respectively, according to the relation^] 



Lad o.5 



[x - 0.5J 
\x + 0.5] 



0.5 
0.5 



(18) 



where [^J an d \%\ denote the integer floor and ceiling 
functions. Additionally, we recall that if S x G [0, 1) is the 
fractional part of x, then [x\ — x — S x . 

Corollary 2 (Bounds on the normalized Nuttall Q): The 
following inequalities can serve as lower and upper bounds 
on the normalized Nuttall Q-function, Qm,n(ch, f3), where 
a,P > and M,N > 0.5, for the case when M > N + 1 
and 5 M = S N (i.e. M-N e W): 



(19) 



with the equalities above being valid only for the case of half- 
odd integer values of M, N. 

Proof: The proof follows immediately from Theorem [2] 

■ 

For the calculation of the bounds in JT9l , the quantities 

QlMJo.s.LJVjo.sK^) and Q[m] .5, Wo.e ( a > P) can be evalu " 
ated exactly by utilizing the results of Corollary Q] Moreover, 
for the case of M, N £ IN, the proposed bounds obtain the 
simplified form 



Q.M,n(0!,P) > QM-0.5,N-0.5(oi, P) 
Q,M,n{01,P) < Q,M+0.5,N+0.5( a >P)- 



(20) 



In Figs. | 2Ja)| and | 2jtb)| the normalized Nuttall Q-fimcu 011 
along with its lower and upper bounds are depicted versus 
P for several values of a and M, respectively, while the 
parameter N is restricted according to the relation N = Al — c 
with c e IN taking values c — 2 in Fig. | 2Ja)| and c = 1, 2, 3 in 
Fig- 1 2jtb)| It is evident, that the bounds proposed in ( fT9l are 
very tight, especially the upper one for Sm(= <5jv) < 0.5 and 
the lower one for 8m(= $n) > 0.5, the latter being the case 
illustrated in Fig. | 2j{b)| 

In order to obtain lower and upper bounds for the standard 
Nuttall Q-function, a similar procedure can be carried out. The 

5 The defining equations of (18) can be easily verified to be valid in any 
arbitrary segment [n, n + 1), where n£Z. 
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Fig. 1. Monotonicity of Qm,n(&i P) with respect to the sum M + N for 
several real values of a , B. 



next theorem will be particularly useful for the fulfillment of 
such a derivation. 

Theorem 3 (Monotonicity of the standard Nuttall Q): The 
standard Nuttall g-f unc ti° n . Q MN (a,/3), where M > 0, 
iV > — 1 and a > 1, p > 0, is strictly increasing with 
respect to the sum M + N, under the requirement of constant 
difference M — N > 1. 

Proof: In Theorem [2] it has been proved that 



d_ 

dv 



Qv+c v=c(a,P) > 



(21) 



where we have substituted v = M + N and c = M — N. 
From (0, d2Tb and after using the quotient rule for partial 
differentiation, we obtain 

^Q^^{a,P)>^ Qs¥ ^{a,P). 



6 
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M=5.0, N=3.0 
M=4.5, Af=2.5 
5.5, N=3.5 
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(a) Q,5,3(a, j3) versus /3 for several values of a. 



1(V 




(b) Qu 2.7(3.5, ft) versus f) for several values of M. 

Fig. 2. Bounds of Qm n(c, 0) for M — N £ IN and several values of 
a, p. 



Since the Nuttall £>-function j s strictly positive, then for a > 1 
it follows that 



e z±£ ^ £ (a,/3)>0 



and the proof is complete. ■ 
Corollary 3 (Bounds on the standard Nuttall Q): The fol- 
lowing inequalities can serve as lower and upper bounds on 
the standard Nuttall (^-function, Q M N (a, B), where a > 1, 
B > and M, N > 0.5, for the case when M > N + 1 and 
Sm = <5 W (i.e. M-N £ M): 



(22) 



with the equalities above being valid only for the case of half- 
odd integer values of M, N. 



Proof: The proof follows immediately from Theorem [3] 

■ 

Similarly to the case of the normalized Nuttall Q, in 
the calculation of the bounds from d22"i l. the quantities 

Q[MW 5 ,iN\ Ja,B) and Q\ M ^,\N^ («. P) can be evalu " 
ated exactly from Theorem Q] Finally, for M, N £ IN, the 
standard Nuttall Q-fimction ca n be simply bounded by 



Qm,n(o:,I3) > QM-0.5,N-oA a > P) 
Qm,n(^,B) < Q M +0.5, N+oA a ^) 



which constitutes the counterpart of 
Nuttall g-fimction. 



(23) 

for the standard 



iii. monotonicity and bounds for the 
Generalized Marcum Q-Function 

Recently, Li and Kam in l46l eq. (11)], following a geo- 
metric approach, presented a novel closed-form formula for 
the evaluation of Q M (a,(3), for the case when M is an odd 
multiple of 0.5 and a > 0, 8 > 0, given by 



Qm((*,P) = ^ erfc 



1 



2q 



V2 

M-1.5 

E 

k=0 
1 



-erfc 



E 



B-a 

~7T 



1 

+ 2 

2 k ^ (k-q)\q\ 



i=0 v ^' 



(-iy 



(24) 



where erfc(z) = (2/y/n) J z e~* dt is the complementary 
error function [3 eq. (7.1.2)]. This representation involves 
only elementary functions and is convenient for evaluation 
both numerically and analytically. For the trivial case when 
a = 0, exact values of the generalized Marcum Q-ftmcuon 
can be obtained from ll46l eq. (12)] 



e M (0,/3)=erfc 



y/2 



_S± M-1.5 



E 



/2tt 



E 

k=0 



2k+l 



-tk-l 



(25) 



Following an algebraic approach, an alternative more compact 
closed-form expression, equivalent to (l24l . can be derived, 
considering the next steps. Particularly, in ||58l eq. (10)] it 
has been proved that the generalized Marcum Q-fimction of 
order m — \i, with m positive integer and < fi < 1 can 
be written in terms of the generalized Marcum g-fimcrion of 
order 1 — // as 



Q m -u(a,B) = e 



m— 1 

E 

71=1 



In-^aB) 



By substituting /i = 0.5 in the above equation and noting that 
for this case the modified Bessel function of the first kind can 
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be replaced by (|7), we obtain 

Qm-0.5 



, 2 

{a,P) = a\l-e 



71=1 



(-2a 2 )" 



(n - fc) n . 



2o 



fc! 



i(2a/3) fe [1 - (-l) fe e 



fc 2a/3] 



(26) 



where once again (m) n denotes the Pochhammer's symbol 
and the term Q 05 (a,f3) can be derived from the definition 
of the generalized Marcum Q-function m <EK by using (3] 
eq. (10.2.14)] as follows 



Go.s(a,/3) 



2 

7T J 



cosh (ax)dx. 



The above integral can be computed in closed-form as 



Go.b(«,j8) 



1 

-erfc 

2 

G(/M 



V2 
a) + GO? 



1 f 
-erfc 

2 

- a) 



/3-a 

~7T 



(27) 



where 2 denotes the Gaussian g-function (or Gaussian 
probability integral) J3] eq. (26.2.3)], defined by Q(z) = 
(1/V2tt) jf z °° e~ t2 l 2 dt. Using ([26) and (27), the generalized 
Marcum Q-function of half-odd integer order can be computed 
for all q > 0, (3 > from the expression 



n-1 



(t» + <3) z 



M-0.5 



(n ^n=L (2a0) k [1 - (-l) fc e 2 ^l 
fc! L J 

+ + g(/3-a), M + 0.5 6M. 



fc=0 



(28) 



We note here that a similar result to ( I28t has been recently 
reported in the literature [48 eq. (16)]. In order to examine the 
special case when a = 0, we first notice that from (@) and (TTSl 
an alternative expression — equivalent to [3, eq. (26.4.25)] — for 
the generalized Marcum 2-function can be derived, written as 



Qm («,/?) 



= P. 2 



^ 2 k k\ 

fc=0 



y 2k r(k + M,^- 

T(k + M) 



a > 0,(3 > 
(29) 



which for integer M falls into the series expansion 
eq. (4)]. Since Q M (a,f3) is a continuous function of a for 
all f3 > and M > 0, the above equation can be extended to 
be asymptotically valid for the case when a — as well, with 
the corresponding limiting value given by 



r 



r(M) 



(30) 



This last result also appears in |[T4l eq. (4.71)], where 
Q M (0,(3) has been derived directly from d3) by applying the 
small argument form of the modified Bessel function. 

It has been proved in |06l that (f24), (f25) along with © 
can define tight upper and lower bounds for the generalized 



Marcum Q-function of integer order. It seems apparent, that 
in order to derive bounds for Q M (a,f3) of real order M, a. 
strict inequality, involving the whole range of M, has to be 
established. Such a generalization concept can be formalized 
through the following theorem. 

Theorem 4 (Monotonicity of the generalized Marcum Q): 
The generalized Marcum g-fu nc ti° n > Q M (a,f3), is strictly 
increasing with respect to its real order M > for all 
a > 0,/3 > 0. 

Proof: Concerning the case when a = 0, we notice that 
can be rewritten as 



Q M (0,/3) = l-P\M,!j 

However, in ll66l eq. (59)] the regularized lower incomplete 
gamma function P(r, x) has been proved to decrease mono- 
tonically with respect to r > for all x > 0. Additionally, for 
a > 0, © implies that the normalized Nuttall Q-function 
with N — M — 1 falls into the generalized Marcum Q- 
function of order M. Nevertheless, according to Theorem [2] 
Qm,m-i(cc, /3) is strictly increasing with respect to 2M — 1 
for M > 0, and the proof is complete. ■ 
The result of Theorem [4] has also recently demonstrated 
by Sun and Baricz in ||62| . where two totally different proofs 
were given. The first one combines the series form of the gen- 
eralized Marcum ^-function presented in (1291 , (130) together 
with the fact that the regularized upper incomplete gamma 
function Q(r,x) = 1 — P(r,x) is strictly increasing with 
respect to r > for each x > 0, originally stated by Tricomi 
in (66 1 . A slightly different analytical proof to this can also be 
found in [60, Th. 1]. The second proof exploits the interesting 
relationship between the generalized Marcum Q-function and 
the reliability function (or CCDF) R of a x 2 random variable 
with 2M DOF and noncentrality parameter a, namely the fact 

that if /3 ~ XiMa 111611 R (P) = Qm(V®i VP)- The interested 
reader is referred to ll62l Th. 3.1] for more information. 

Recalling the relation between the normalized Nuttall and 
the generalized Marcum Q-functions, that is (@), Fig. [l|fa)| 



verifies graphically the results of Theorem |4] since it actually 
depicts Q M (a, f3) versus the term 2M — 1. 

Corollary 4 (Bounds on the generalized Marcum Q): The 
following inequalities can serve as lower and upper bounds 
on the generalized Marcum 2-f unc ri on Q M (a,(3) of real 
order M > 0.5 for all a > 0, /3 > 0. 



(31) 



with the equalities above being valid only for the case of half- 
odd integer values of M. 

Proof: The proof follows immediately from Theorem [4] 

■ 

In CorollaryH] the quantities QimJo b ( a ' @) an< ^ 6[mi . 5 ( a > P) 
can be evaluated exactly either from (T24l . (f25) or ( |27] |. d30l) . 
while for Mgl d3l1 > reduces to 

QM-o.5( a iP) < GmK/ 3 ) < Qm+oA^P) 

which comes as a complement to the inequalities of (|20l 
and (|23l . This last result was originally demonstrated in 1461 
eq. (16)], where the authors following a geometric approach 
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(a) Qi(ct, p) versus /3 for several values of a. 




(b) Qm (2-5, 0) versus j3 for several values of M. 

Fig. 3. Bounds of Q M (a, 0) for seve ral real values of a, /3 and M. (In |(a)| 
bounds proposed by Li and Kara, 1461 ). 



proposed tight lower and upper bounds for the generalized 
Marcum 2-f unct i° n of integer order M, which have been 
proved to outperform other existing ones. This can be easily 
verified from Fig.|3 |Ia)| where Q i (a 1 0) has been plotted versus 
j3 for several values of a. Therefore, for the case of real M, 
one can expect even further enhancement in the strictness 
of either the lower bound (for Sm > 0.5) or the upper one 
(for Sm < 0.5). This is clearly depicted in Fig. 1 3fb)| where 
the curves Q 2 5 (2.5, /3) and Q 8 5 (2.5, (3) constitute very tight 
lower and upper bounds of Q 2 ^(2.5, (3) and <2 8 3 (2.5, /3), 
respectively, for all range of (3. 

IV. Conclusion 

Applicable monotonicity criteria were established for the 
normalized and standard Nuttall and the generalized Marcum 



Q-functions. Specifically, it was proved that the two Nuttall Q- 
functions are strictly increasing with respect to the real sum 
M + N for the case when M > N + 1, while the generalized 
Marcum (^-function increases monotonically with respect to 
its real order M. Additionally, novel closed-form expressions 
for both types of the Nuttall (^-function we re given for the case 
when M, N are odd multiples of 0.5 and M > N. Regarding 
the generalized Marcum Q-fimction G f half-odd integer order, 
an alternative more compact closed-form expression, equiv- 
alent to the already existing one, was derived. By exploiting 
these results, novel lower and upper bounds were proposed for 
the Nuttall Q-functions when M > N + 1, while the recently 
proposed bounds for the generalized Marcum Q-f\mcl\on of 
integer M, were appropriately utilized in order to extend their 
validity over real values of M. 

References 

[1] J. I. Marcum, "A statistical theory of target detection by pulsed radar," 
IRE Trans. Inf. Theory, vol. IT-6, no. 2, pp. 59-267, Apr. 1960, 
(Originally published by the Rand Corporation, Santa Monica, CA, 
Research Memorandum RM-754, Dec. 1947, and RM-753, Jul. 1948). 

[2] H. Nuttall, "Some integrals involving the Q-function," Naval Underwater 
Systems Center, New London, CT, Tech. Rep. 4297, Apr. 1972. 

[3] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions 
with Formulas, Graphs, and Mathematical Tables, 9th ed. New York, 
NY: Dover Publications, Inc., 1972. 

[4] M. K. Simon, "The Nuttall Q function — its relation to the Marcum 
Q function and its application in digital communication performance 
evaluation," IEEE Trans. Commun., vol. 50, no. 11, pp. 1712-1715, 
Nov. 2002. 

[5] H.-C. Yang and M.-S. Alouini, "Closed-form formulas for the outage 
probability of wireless communication systems with a minimum signal 
power constraint," IEEE Trans. Veh. Technol., vol. 51, no. 6, pp. 1689— 
1698, Nov. 2002. 

[6] L. Yang and M.-S. Alouini, "On the average outage rate and aver- 
age outage duration of wireless communication systems with multiple 
cochannel interferers," IEEE Trans. Wireless Commun., vol. 3, no. 4, pp. 
1142-1153, Jul. 2004. 
[7] J. M. Romero-Jerez and A. J. Goldsmith, "Antenna array processing 
in fading and interference: An interference-cancellation vs. diversity 
comparative performance," in Proc. IEEE Int. Conf. on Communications 
(ICC'07), Glasgow, Scotland, Jun. 24-28, 2007, pp. 1108-1113. 
[8] M. Kang and M.-S. Alouini, "Largest eigenvalue of complex Wishart 
matrices and performance analysis of MIMO MRC systems," IEEE J. 
Sel. Areas Commun., vol. 21, no. 3, pp. 418^126, Apr. 2003. 
[9] S. Jin, M. R. Mckay, X. Gao, and I. B. Collings, "MIMO multichannel 
beamforming: SER and outage using new eigenvalue distributions of 
complex noncentral Wishart matrices," IEEE Trans. Commun., vol. 56, 
no. 3, pp. 424-^134, Mar. 2008. 

[10] A. Maaref and S. Aissa, "Joint and marginal eigenvalue distributions 
of (non)central complex Wishart matrices and PDF-based approach for 
characterizing the capacity statistics of MIMO Ricean and Rayleigh 
fading channels," IEEE Trans. Wireless Commun., vol. 6, no. 10, pp. 
3607-3619, Oct. 2007. 

[11] K.-P. Ho and H.-C. Wang, "Soft-decoding of low-density parity 
check codes for DPSK signals," in Proc. Optical Fiber Communica- 
tion/National Fiber Optic Engineers Conf. (OFC/NFOEC'07), Anaheim, 
CA, Mar. 25-29, 2007, pp. 1-3. 

[12] J. I. Marcum, "Table of Q functions," Rand Corporation, Santa Monica, 
CA, Research Memorandum RM-339, Jan. 1950. 

[13] A. H. Nuttall, "Some integrals involving the Q M function," IEEE Trans. 
Inf. Theory, vol. 21, no. 1, pp. 95-96, Jan. 1975, (Originally published 
by the Naval Underwater Systems Center, New London, CT, Tech. Rep. 
4755, May 1974). 

[14] M. K. Simon and M.-S. Alouini, Digital Communication over Fading 
Channels, 2nd ed. Hoboken, NJ: Wiley-Interscience, 2005. 

[15] D. A. Shnidman, "The calculation of the probability of detection and 
the generalized Marcum (J-function," IEEE Trans. Inf. Theory, vol. 35, 
no. 2, pp. 389^100, Mar. 1989. 



KAPINAS et at: ON THE MONOTONICITY OF THE GENERALIZED MARCUM AND NUTTALL G-FUNCTIONS 



[16] P. Swerling, "Probability of detection for fluctuating targets," IRE Trans. 

Inf. Theory, vol. IT-6, no. 2, pp. 269-308, Apr. 1960, (Originally 

published by the Rand Corporation, Santa Monica, CA, Research 

Memorandum RM-1217, Mar. 1954). 
[17] C. W. Helstrom, Statistical Theory of Signal Detection, 2nd ed. New 

York, NY: Pergamon Press, 1968. 
[18] H. L. V. Trees, Detection, Estimation, and Modulation Theory, Part I: 

Detection, Estimation and Linear Modulation Theory. New York, NY: 

Wiley-Interscience, 2001. 
[19] M. K. Simon and M.-S. Alouini, "A unified approach to the probability 

of error for noncoherent and differentially coherent modulations over 

generalized fading channels," IEEE Trans. Commun., vol. 46, no. 12, 

pp. 1625-1638, Dec. 1998. 
[20] , "A unified performance analysis of digital communication 

with dual selective combining diversity over correlated Rayleigh and 

Nakagami-m fading channels," IEEE Trans. Commun., vol. 47, no. 1, 

pp. 33^13, Jan. 1999. 
[21] , "Some new results for integrals involving the generalized Marcum 

Q function and their application to performance evaluation over fading 

channels," IEEE Trans. Wireless Commun., vol. 2, no. 4, pp. 611-615, 

Jul. 2003. 

[22] S. Gaur and A. Annamalai, "Some integrals involving the 

Q m (a^/x, by/x) with application to error probability analysis of 

diversity receivers," IEEE Trans. Veh. Technol., vol. 52, no. 6, pp. 

1568-1575, Nov. 2003. 
[23] A. Annamalai, C. Tellambura, and V. K. Bhargava, "A general method 

for calculating error probabilities over fading channels," IEEE Trans. 

Commun., vol. 53, no. 5, pp. 841-852, May 2005. 
[24] F. F. Digham, M.-S. Alouini, and M. K. Simon, "On the energy detection 

of unknown signals over fading channels," IEEE Trans. Commun., 

vol. 55, no. 1, pp. 21-24, Jan. 2007. 
[25] A. Pandharipande and J. -P. M. G. Linnartz, "Performance analysis of 

primary user detection in a multiple antenna cognitive radio," in Proc. 

IEEE Int. Conf. on Communications (ICC'07), Glasgow, Scotland, Jun. 

24-28, 2007, pp. 6482-6486. 
[26] Y. Zhu, P. Y. Kam, and Y. Xin, "On the mutual information distribution 

of MIMO Rician fading channels," IEEE Trans. Commun., to be 

published. 

[27] M. K. Simon, Probability Distributions Involving Gaussian Random 
Variables: A Handbook for Engineers and Scientists. New York, NY: 
Springer Science-Business Media, LLC, 2006. 

[28] D. A. Shnidman, "Two noncentral chi-square generalizations," IEEE 
Trans. Inf. Theory, vol. 42, no. 1, pp. 283-285, Jan. 1996. 

[29] J. G. Proakis, Digital Communications, 4th ed. New York, NY: 
McGraw-Hill, 2000. 

[30] J. Omura and T. Kailath, "Some useful probability distributions," Sys- 
tems Theory Lab., Stanford Univ., Stanford, CA, Tech. Rep. 7050-6, 
Sep. 1965. 

[31] M. Schwartz, W. R. Bennett, and S. Stein, Communication Systems and 

Techniques. New York, NY: McGraw-Hill, 1966. 
[32] M. K. Simon and M.-S. Alouini, "A simple single integral representation 

of the bivariate Rayleigh distribution," IEEE Commun. Lett., vol. 2, no. 5, 

pp. 128-130, May 1998. 
[33] G. V. Weinberg, "Poisson representation and Monte Carlo estimation of 

generalized Marcum Q-function," IEEE Trans. Aerosp. Electron. Syst., 

vol. 42, no. 4, pp. 1520-1531, Oct. 2006. 
[34] G. H. Robertson, "Computation of the noncentral chi-square distribu- 
tion," Bell Syst. Tech. J., vol. 48, pp. 201-207, Jan. 1969. 
[35] G. M. Dillard, "Recursive computation of the generalized Q function," 

IEEE Trans. Aerosp. Electron. Syst., vol. AES-9, no. 4, pp. 614-615, 

Jul. 1973. 

[36] D. A. Shnidman, "Efficient evaluation of probabilities of detection and 

the generalized Q-function," IEEE Trans. Inf. Theory, vol. IT-22, no. 6, 

pp. 746-751, Nov. 1976. 
[37] S. S. Rappaport, "Computing approximations for the generalized Q 

function and its complement," IEEE Trans. Inf. Theory, vol. IT-17, no. 4, 

pp. 497^98, Jul. 1971. 
[38] C. W. Helstrom, "Computing the generalized Marcum Q-function," IEEE 

Trans. Inf. Theory, vol. 38, no. 4, pp. 1422-1428, Jul. 1992. 
[39] G. V. Weinberg and L. Panton, "Numerical estimation of Marcum's Q- 

function using Monte Carlo approximation schemes," Defence Science 

and Technology Organisation, Edinburgh, Australia, Research Rep. 

DSTO-RR-0311, Apr. 2006. 
[40] N. Ding and H. Zhang, "A flexible method to approximate Marcum 

Q-function based on geometric way of thinking," in Proc. Int. Symp. 

on Communications, Control and Signal Processing (ISCCSP'08), St. 

Julians, Malta, Mar. 12-14, 2008, pp. 1351-1356. 



[41] D. E. Johansen, "New techniques for machine computation of the Q- 
functions, truncated normal deviates, and matrix eigenvalues," Applied 
Research Laboratory, Sylvania Electronic Systems, Waltham, MA, Tech. 
Rep. AFCRL-556, Jul. 1961. 

[42] S. Pari, "A new method of calculating the generalized Q function," IEEE 
Trans. Inf. Theory, vol. IT-26, no. 1, pp. 121-124, Jan. 1980. 

[43] P. E. Cantrell, "On the calculation of the generalized Q-function via 
Part's method," IEEE Trans. Inf. Theory, vol. IT-32, no. 6, pp. 817-824, 
Nov. 1986. 

[44] M. K. Simon, "A new twist on the Marcum Q-function and its applica- 
tion," IEEE Commun. Lett, vol. 2, no. 2, pp. 39-41, Feb. 1998. 

[45] A. H. M. Ross, "Algorithm for calculating the noncentral chi-square 
distribution," IEEE Trans. Inf. Theory, vol. 45, no. 4, pp. 1327-1333, 
May 1999. 

[46] R. Li and P. Y. Kam, "Computing and bounding the generalized Marcum 
Q-function via a geometric approach," in Proc. IEEE Int. Symp. on 
Information Theory (ISIT'06), Seattle, WA, Jul. 9-14, 2006, pp. 1090- 
1094. 

[47] G. T. F. de Abreu, "Jensen-Cotes upper and lower bounds on the 
Gaussian Q-function and related functions," IEEE Trans. Commun., 
submitted for publication. 

[48] Y. Sun and S. Zhou, "Tight bounds of the generalized Marcum Q- 
function based on log-concavity," in Proc. IEEE Global Telecommu- 
nications Conf. (GLOBECOM'08), New Orleans, LA, Nov. 30-Dec. 4, 
2008, pp. 1-5. 

[49] G. E. Corazza and G. Ferrari, "New bounds for the Marcum g-function," 

IEEE Trans. Inf. Theory, vol. 48, no. 11, pp. 3003-3008, Nov. 2002. 
[50] P. Y. Kam and R. Li, "Computing and bounding the first-order Marcum 

Q-function: A geometric approach," IEEE Trans. Commun., vol. 56, 

no. 7, pp. 1101-1110, Jul. 2008. 
[51] X. Zhao, D. Gong, and Y. Li, "Tight geometric bound for Marcum Q- 

function," Electron. Lett, vol. 44, no. 5, pp. 340-341, Feb. 2008. 
[52] R. Li and P. Y. Kam, "Generic exponential bounds on the generalized 

Marcum Q-function via the geometric approach," in Proc. IEEE Global 

Telecommunications Conf. (GLOBECOM'07), Washington, DC, Nov. 

26-30, 2007, pp. 1754-1758. 
[53] M. K. Simon and M.-S. Alouini, "Exponential-type bounds on the 

generalized Marcum Q-function with application to error probability 

analysis over fading channels," IEEE Trans. Commun., vol. 48, no. 3, 

pp. 359-366, Mar. 2000. 
[54] A. Annamalai and C. Tellambura, "Cauchy-Schwarz bound on the 

generalized Marcum Q-function with applications," Wirel. Commun. 

Mob. Comput., vol. 1, no. 2, pp. 243-253, Mar. 2001. 
[55] R. F. Mills and G. E. Prescott, "A comparison of various radiometer 

detection models," IEEE Trans. Aerosp. Electron. Syst., vol. 32, no. 1, 

pp. 467^73, Jan. 1996. 
[56] C. W. Helstrom, "Approximate inversion of Marcum's Q-function," IEEE 

Trans. Aerosp. Electron. Syst., vol. 34, no. 1, pp. 317-319, Jan. 1998. 
[57] S. P. Herath and N. Rajatheva, "Analysis of equal gain combining 

in energy detection for cognitive radio over Nakagami channels," in 

Proc. IEEE Global Telecommunications Conf. (GLOBECOM'08), New 

Orleans, LA, Nov. 30-Dec. 4, 2008, pp. 1-5. 
[58] G. Di Blasio and A. Neri, "Modified Q-functions and their use in 

detection analysis," IEEE Trans. Inf. Theory, vol. 37, no. 4, pp. 1123- 

1142, Jul. 1991. 

[59] P. B. Patnaik, "The non-central x 2 ~ an d -F-distributions and their 
applications," Biometrika, vol. 36, no. 1/2, pp. 202-232, Jun. 1949. 

[60] B. K. Ghosh, "Some monotonicity theorems for x 2 , F and t distri- 
butions with applications," /. Roy. Statist. Soc. B, vol. 35, no. 3, pp. 
480-492, 1973. 

[61] N. M. Temme, "Asymptotic and numerical aspects of the noncentral 
chi-square distribution," Comput. Math. Appi, vol. 25, no. 5, pp. 55-63, 
1993. 

[62] Y. Sun and A. Baricz, "Inequalities for the generalized Marcum Q- 
function," Appl. Math. Comput, vol. 203, no. 1, pp. 134-141, Sep. 2008. 

[63] A. Baricz and Y. Sun, "New bounds for the generalized Marcum Q- 
function," IEEE Trans. Inf. Theory, vol. 55, no. 7, pp. 3091-3100, Jul. 
2009. 

[64] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and 

Products, 6th ed. San Diego, CA: Academic Press, 2000. 
[65] A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev, Integrals and 

Series, Vol. 1: Elementary Functions, 3rd ed. New York, NY: Gordon 

and Breach Science Publications, 1992. 
[66] F. G. Tricomi, "Sulla funzione gamma incompleta," Ann. Mat. Pura 

Appl, vol. 31, no. 1, pp. 263-279, Dec. 1950. 



10 



PUBLISHED IN IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 55, NO. 8, AUGUST 2009 



[67] F. Qi, "Monotonicity results and inequalities for the gamma and incom- 
plete gamma functions," Math. Inequal. Appl, vol. 5, no. 1, pp. 61-68, 
2002. 

[68] E. C. Titchmarsh, The Theory of Functions, 2nd ed. London, UK: 
Oxford Univ. Press, 1939. 



Vasilios M. Kapinas (S'07-M'09) was born in Thessaloniki, Greece, in May 
1976. He received the diploma degree in electrical and computer engineering 
from Aristotle University of Thessaloniki, Greece, in 2000. Since 2005, he has 
been working toward the Ph.D. degree in telecommunications engineering. 

His current research interests include wireless communication theory and 
digital communications over fading channels, giving special focus to space- 
time block coding techniques. 



Sotirios K. Mihos was born in Thessaloniki, Greece, in April 1984. He is 
an undergraduate student at the Aristotle University of Thessaloniki, Greece, 
where he is working toward the diploma degree in electrical and computer 
engineering. 

His research interests span a wide range of subject areas including 
computer science, electronics and automatic control, with a special focus on 
their relationship to pure mathematics. 



George K. Karagiannidis (M'97-SM'04) was born in Pithagorion, Samos 
Island, Greece. He received the University and Ph.D. degrees in electrical 
engineering from the University of Patras, Patras, Greece, in 1987 and 1999, 
respectively. From 2000 to 2004, he was a Senior Researcher at the Institute 
for Space Applications and Remote Sensing, National Observatory of Athens, 
Greece. In June 2004, he joined Aristotle University of Thessaloniki, Thessa- 
loniki, Greece, where he is currently an Assistant Professor in the Electrical 
and Computer Engineering Department. His current research interests include 
wireless communication theory, digital communications over fading channels, 
cooperative diversity systems, cognitive radio, satellite communications, and 
wireless optical communications. 

He is the author or coauthor of more than 80 technical papers published 
in scientific journals and presented at international conferences. He is also 
a coauthor of two chapters in books and a coauthor of the Greek edition of 
a book on mobile communications. He serves on the editorial board of the 
EURASIP Journal on Wireless Communications and Network- 
ing. 

Dr. Karagiannidis has been a member of Technical Program Committees 
for several IEEE conferences. He is a member of the editorial boards of the 
IEEE Transactions on Communications and the IEEE Communica- 
tions LETTERS. He is co-recipient of the Best Paper Award of the Wireless 
Communications Symposium (WCS) in IEEE International Conference on 
Communications (ICC 07), Glasgow, U.K., June 2007. He is a full member 
of Sigma Xi. 



